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Temperature  Distribution  in  the  Boundary 
Layer  of  an  Airplane  Wing  with  a Line  Source  of 
Heat  at  the  Stagnation  Edge  ~ 

Symmetric  Wing  in  Symmetric  Plow 

by 

Chia-Shun  Yih,  Jack  E.  Cermak,  and  Richard  T.  Shen 

Abstract 

Given  a symmetric  airplane  wing  in  a symmetric  flow 
with  an  insulated  surface  and  a line  source  of  heat  at  the 
stagnation  edge,  it  is  proposed  to  calculate  the  temperature 
distribution  in  the  boundary  layer,  the  minor  effect  of  free 
convection  being  neglected.  With  the  velocity  distribution 
given  by  previous  writers  in  the  form  of  a power  series  of 
the  curvilinear  abscissa,  a similar  expansion  is  assumed  for 
the  temperature  distribution  and  the  resulting  ordinary 
differential  equations  for  the  functions  occurring  in  the 
coefficients  of  the  powers  are  solved  numerically.  Since 
these  functions  do  not  depend  on  the  form  of  the  symmetric 
wing  section,  they  are  universal  functions  that  can  be 
applied  to  symmetric  wing  sections  of  any  form  in  flow 
without  yaw,  for  calculating  the  temperature  distribution 
in  the  boundary  layer,  A sample  calculation  of  the  tempera- 
ture distribution  within  the  boundary  layer  of  a circular 
cylinder  in  symmetric  flow  is  presented.  The  results  will 
also  be  found  to  be  useful  when  an  unsyusnetric  wing  or  a 
syrmetric  wing  with  yaw  is  considered. 


1.  Introduction 


As  a first  seep  toward  tho  investigation  of  the  possi- 
bility of  preventing  icing  on  airplane  wings  by  a line  source 
of  heat  and  ma3s,  one  considers  a 3ynmetric  wing  in  synr.etric 
flow  with  an  insulated  surface  and  a line  source  of  heat 
(heating  filament)  at  the  stagnation  edge,  and  seeks  to 
find  the  temperature  distribution  in  the  boundary  layer. 

If  the  effects  of  free  convection  and  compressibility  are 
neglected,  the  velocity  distribution  in  the  boundary  layer 
is  already  given  by  previous  writers  in  the  form  of  a 
po'fer  series  of  the  curvilinear  abscissa  — Goldstein 
(1:151 )•  A similar  series  can  be  assumed  for  the  temper- 
ature distribution.  When  this  series  is  substituted  in 
the  boundary- layer  equation  of  diffusion,  ord?r»ry  differ- 
ential equations  for  the  functions  occurring  in  the  co- 
efficients of  the  powers  will  be  obtained.  These 
equations,  which  do  not  depend  on  the  form  of  the  wing 
section,  will  be  solved  numerically.  The  solutions  provide 
universal  functions  which  are  not  only  sufficient  for 
calculating  the  temperature  distribution  in  the  boundary 
layer  of  a symmetric  wing  of  any  form  in  symmetric  flow, 
but  also  useful  when  an  unsymmetric  wing  or  a symmetric 
wing  with  yaw  is  considered.  Calculations  for  cases  of 
unsymmetric  flow  are  not  included  in  this  report# 

2.  Velocity  Distribution  in  ths  Boundary  Layer 

In  a plane  perpendicular  to  the  length  of  the  wing, 
the  trace  of  the  stagnation  edge  will  be  taken  to  be  the 
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origin,  r will  be  measured  along  the  boundary  of  the 
wing,  and  y^  will  be  measured  in  a direction  normal  tc 
that  or  x^  (see  Fig.  1 )#  Denoting  by  u-^  and 
the  velocity  components  in  the  directions  of  x^  and  y 
respectively,  by  UQ  the  velocity  of  approach,  and  by- 
IT^  the  potential  velocity  outside  of  the  boundary  layer 
in  the  direction  of  x^  , one  has  the  following  boundary- 
layer  equation  of  motion: 


u,|^+  v.^  = u + ' 

‘3x-  '3y.  U'dx,  Vay* 


X, 

Ux  being  a function  of  x^  only,  and  y being  the 
kinematic  viscosity.  In  order  to  eliminate  y explicitly 
and  to  deal  only  with  dimensionless  quantities,  one  makes 
the  following  substitutions: 


(1) 


u.  . u' 


x-D-y-^y,^  »•■*>  v ’ u-u, 

where  D is  a reference  length  of  the  wing  section,  and 
transforms  (1)  into 

uu,  + vuv  - U U/+  u 


(2) 


'yy 


(3) 


v;here  \j'~  — — . 

dx 


The  equation  of  continuity  can  be  written  in  the 
dimensionless  form 

ux  * Vy  *0  (4) 

which  permits  the  use  of  a dimensionless  stream- function 
such  that  bt  = vpy  , V = -<px 

iT  For  the  meaning  of  symbols  see  LIST  OP  SYMBOLS  on  page  11. 


Thus  (3) 


becomes 


-v>„4>yy  = UU'+<f<yyy 

l/ith  the  dimensionless  potential  velocity  given  by 


U = <ajX  q3  x5  + asx5+ 

where  at , C^,  d5»  are  dimensionless  quantities  depending 
only  on  the  form  of  the  wing  section,  and  where  U is 
an  odd  function  of  x due  to  the  definitions  of  U and 
x and  due  to  symmetry,  one  may  assume 

$ = </j  * 4 <^X3  + <^X5  + ^x7  + 


where 


* Va,,Yj  vs;  * 


where  f , g , h , etc.  are  functions  of  the  new  variable 

n =y  Ja, 

Prom  (7)  and  (8)  one  has 

g'y+-^h;  ) xr+ 

u , 


- - ^9,+  ■ * 4 ) 

A series  of  equations  in  f , g , h etc.  are  furnished 

by  (5)  whose  boundary  conditions  are 

Uo)  =9<0)  =hCO)  - = 0 (for  oil  subscripts; 


fJ/(oo)  * etc. 

The  solutions  for  f^  , f^  , g^  and  h^  are  given  by 
Goldstein  (1:151)  and  Schlichting  (2:122). 


(5) 

(6) 

(7) 

W 

(9) 

(10) 

(ID 

(12) 

(13) 


3»  Temperature  Distribution  in  the  Boundary  Layer 


As  the  fluid  passes  the  heat  source,  it  will  carry 
away  some  heat  and  cause  the  tenperature  in  the  thermal 
boundary  layer  to  rise,  while  the  temperature  outside 
of  the  layer  remains  practically  the  same  as  the  ambient 
temperature  TQ,  The  amount  of  heat  carried  into  the 
boundary  layer  is  a measure  of  the  strength  of  the  line 
source.  Denoting  this  strength  by  2H  , and  the  temper- 
ature at  any  point  by  T , one  has 

H =pcp  [ u(  (T-T0)dy, 
o 

where  p is  the  density  and  c^  is  the  specific  ho.Qt  at 
constant  pressure  of  the  fluid.  If  the  wing  has  an  insu- 
lated surface,  H should  be  a constant  independent  of 
x^.  Using  the  dimensionless  parameters 

e.Irli  , u=-£- 

T„  U„ 

and  the  new  variable  given  in  (9),  one  can  write  (14)  as 


(14) 


US) 


or 


# I ^ 

h \fo{  r°* 

Neglecting  the  heat  generated  by  compressibility  and 
by  viscous  shear,  the  boundary-layer  equation  of  heat 
diffusion  can  be  wiittcr. 


(16) 


1(6,+  V0y 


_L  G 

<r 


(17) 


where  o'  - }L-  is  the  Prandtl  number,  o C being  the  thermal 


diffuaivity.  Remembering  (10)  and  (16),  one  assumes 


0 - 
where 


H r— - 

/0CpTeuo  d ^ (0o+  e*x44 ' 


(18) 


2^ 

e„-  F.(r|),  e,=  p*(n) » [<Vi) + Jt"*'  1^r<  191 

i I o 

and  R = , Substituting  (18)  into  (17),  and  equating 

the  coefficients  of  equal  powers  of  x on  both  aider,  one 
has  the  following  series  of  ordinary  differential  equations: 


iFc'/”-(f,/F„+f,F.') 


( r’o 
I r» 


zh  = f,  f;-  ( V F.  + 3fj F.' ) 

9 ^ -3f '6*-  F.  + 59sF;  ) 

j:Hf-jf'H4-fH4-(h;F0+shsF,'-  8f3F2') 


(20) 

(21) 

(22) 

(23) 


where  the  primes  denote  differentiation  with  respect  to 
• The  boundary  conditions  at  the  insulated  surface 
of  the  wing  (where  I'j-Q  ) are 

F.'fO)  = F2'(0)  - G^COI  = h4(0)' =0  (21t) 

while  those  at  points  outside  of  the  thermal  boundary 
layer  can  oe  written 

F (oc)  ~ =■  * » 0 * 2 ^ 

The  integral  condition  imposed  by  (16)  can  be 
decomposed  after  a straight-forward  calculation  with 
(10)  and  (l8),  into  the  following  integral  conditions: 


Jw'<*0- 

c 

l ( ^ ) d *9  * 


(26) 

(27) 


(28) 


(29) 


j/ls'li+f,'G4)dn=o 

Cf» 

/ (3f,/H4+3h5/F.  + 8fJ/F1)rf^  -O 

It  will  now  be  shown  that  by  virtue  of  (21)  to  (25)  the 
conditions  (2?)  to  (29)  are  always  satisfied.  For 
instance# by  virtue  of  (21)  one  can  write 

2 ( flte + C F,-)  = i>F2"+ + C f;  ) + 3<:  5 f0  + f3  f;  ; 

integration  of  which  yields 

+ X,  f.  17- 0 

by  (24),  (25)  and  (12),  so  that  (27)  is  satisfied.  If 
in  virtue  of  (22)  f,nd  (23),  one  writes 

4 (9jF.  +f,G4)  =£G"+0.g1H'Ga'>  +5’(<35f«,+  gs'FJ 
3 (Jf‘  4)+3hsF.+  SF/F^'j  F/h^I  + s(hsF.+h'sF. )+8 Cf,  F/+f'F) 

(28)  and  (29)  can  be  similarly  demonstrated  upon 
integration.  In  fact,  all  the  integral  conditions 
subsequent  to  (29)  are  always  satisfied  if  the  differ- 
ential equations  and  the  boundary  conditions  (24)  and 
(25)  are  satisfied.  Thus,  only  (26)  remains,  which 
takes  the  place  of  (16). 

Integrating  (20),  one  obtains 

— P'=  - f P 
cr  ° , c 
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(30) 


S\ 


the  constant  of  integration  being  zero  since  Fc  (0)  = 
f^O)  = 0«  A second  integration  gives 

F = k 


(3D 


where  K is  determined  from  (26)  to  be 

0 

Since  the  functions  fx  , f^  , , h^  as  well  as  their 

derivatives  are  tabulated  by  Goldstein  (1:151)  and 
Schlichting  (2:122),  one  can  calculate  Fq  numerically. 

The  results  are  given  in  Table  1 for  CT  ~ 0.73* 

In  solving  (21)  numerically  one  takes  a trial  value 
for  (0)  * and  utilizing  the  boundary  condition 

(0)  =0  , proceeds  to  calculate  the  function  F2*  If 
the  condition  at  infinity  is  not  satisfied,  another  trial 
value  for  F^  (°)  is  taken,  until  the  condition  F2  ( oo  ) = o 
Is  satisfied.  This  method  of  solution  also  applies  to  (22) 
and  (23),  Values  for  g^  and  h^  and  their  fir3t  two 
derivatives  for  intermediate  values  of  fy  not  entered  in 
Table  III  of  Schlichting  (2:122)  were  obtained  by  inter- 
polation according  to  the  method  given  by  Sokolnikoff 
(3:55D.  The  results  for  F2,  , and  are  given 

respectively  in  Tables  2,  3>  and  k.  The  functions  F0  , 

F2  * &nC^  Hk  are  rePresentec^  graphically  in  Figure  2, 

k*  A Particular  Example 

In  order  to  illustrate  the  application  of  the  fore- 


going results,  temperature  distributions  within  the 


boundary  layer  of  a right  circular  cylinder  are  calculated. 
With  the  longitudinal  axis  of  the  cylinder  oriented 
perpendicular  to  the  direction  of  the  undisturbed  velocity 
field,  radial  temperature  distributions  are  computed  for 
each  10°  of  central  angle  — the  angle  being  measured 
from  the  stagnation  line  — from  10°  to  and  including  80°. 


The  dimensionless  quantity 


9pc  pT0U.  D 
HvTr 


is 


calculated  as  a function  of  for  each  of  the  values  of 

x chosen  according  to  the  foregoing  paragraph.  By  (18), 

A is  found  to  be 


£ “ 7r==^(  + ) (33) 

i/a,  x 2 0 2 * 

and  is  computed  for  various  values  of  rj  at  a particular 
value  of  x by  use  of  (19)  and  Tables  1 through  4*  The 
constants  and  0^  are  taken  to  be  those 

experimentally  determined  by  Kiemenz  — see  Soldstein 
(1:150)  --  under  the  following  conditions: 


U0  = 19 .2  cm/sec 

D = 9.74  cm 
P.  = 1.85  x 10^; 

therefore,  the  values  for  ^3  given  in  Table  5 are 
applicable  for  values  of  R near  1.85  x 10^. 

Figure  3 shows  graphically  how  p variee  with  Kj  for 
tne  various  values  of  x . Having  obtained  values  for  fi  , 
the  temperature  distribution  may  be  calculated  for 
particular  ambient  conditions  and  fluid  properties,  heat 
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3ource  strength,  and  cylinder  diameter  by  means  of  (15) 
and  (16),  for  K ~ 1.85  X tO* 

Corresponding  calculations  for  a symmetrical  airplane 
wing  section  in  symmetric  flow  may  be  carried  out  after  the 
constants  Qf  , and  CX^  have  been  determined  for  the 

wing. 
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K 
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Uo 
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h'  h4 


reference  length  of  the  wing  section 
functions  of  rj 

one-half  the  strength  of  heat  source 
constant 

Reynolds  number  D 

V 

temperature  at  any  point 
ambient  fluid  temperatui  e 
dimensionless  velocitv  parameter 

Uo 

velocity  of  approach 

potential  velocity  outside  of  the  boundary 
layer  in  the  direction  of  x^ 

dimensionless  constants  for  a particular 
shape  of  wing 


c 

p 

specific  heat  of  the  fluid  at  constant 
pressure 

f ^ t htj 

functions  of 

u 

dimensionless  velocity  parameter 

U1 

velocity  component  in  the  direction  of 

X1 

V 

VI 

dimensionless  velocity  parameter  1/  P—  1/ 

Vucv  ' 

velocity  component  in  the  direction  of 

yl 

X 

dimensionless  length  parameter.^! 

*1 

distance  measured  along  the  boundary  of 
wing 
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y 

dimensionless  length  parameter  y 

’ VD  ' 

distance  measured  normal  to  the  direction  of 
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thermal  diffuaivity 
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dimensionless  temperature  parameter 
boundary  layer  thickness 
dimensionless  length  parameter 
dimensionless  tenperature  parameter 
kinematic  viscosity  of  the  fluid 
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dimensionless  strecra-function 
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Table  2 --  Values  of  F2,  F 2 * and  F2. 
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-0.0559 

1.0 
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0.0221 
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0.0011 

-0.0030 

2.3 

0.1226 

0.0916 

-0.1578 

5.8 
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Table  4 

--  Values  of  Hj^, 
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0.0779 
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0.0798 
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0.0318 

0.0262 

0.0212 
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0.0034 
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4.378 

2.200 

1.475 
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1.116 
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0.675 
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0.889 
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2.017 
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0.651 
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0.594 

1.1 

3.750 
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0.589 
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1.6 
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0.562 
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0.629 
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0.524 
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0.555 
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2.0 
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1.839 
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0.479 

0.454 

0.453 
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1.649 
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0.392 
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0.362 
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0.300 
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0.859 

0.454 
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0.738 
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0.288 

0.248 
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0.245 
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2.9 
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0.335 
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0.219 
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0.164 
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3.1 
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0.241 
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0 , 22b 
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0.117 

0.119 

0.130 

0.150 

0.177 

3.4 

0.254 

0.138 

0.107 

0.099 

0.101 

0.112 

0.130 

0.155 
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Figure  3 --  Graphs  of  (9  for  a right  circular  cylinder  in  Byaaaotric  flow. 


